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ABSTRACT 

Based on a suite of Monte Carlo simulations, I show that a stellar-mass- 
dependent lifetime of the gas disks from which planets form can explain the 
lack of hot Jupiters/close-in giant planets around high-mass stars and other key 
features of the observed semimajor axis distribution of radial velocity-detected 
giant planets. Using reasonable parameters for the Type II migration rate, re- 
gions of planet formation, and timescales for gas giant core formation, I construct 
synthetic distributions of Jovian planets. A planet formation/migration model 
assuming a stellar mass-dependent gas disk lifetime reproduces key features in 
the observed distribution by preferentially stranding planets around high-mass 
stars at large semimajor axes. 

Subject headings: Stars: planetary systems: formation 



Introduction 



The semimajor axis distribution of Jovian-mass planets discovered in radial velocity 
surveys reveals striking trends (Figured! top panels). Many Jovian planets around solar-mass 
stars have semimajor axes a p > 0.5 AU. Many also hav e & p < 0.1-0.2 AU ('h ot Jupiters' 



and few have intermediate values ('the period valley', ICumming et al.ll2008l ). While hot 
Jupiters comprise ~ 20% of planets around < 1.5 M stars, surveys have yet to detect hot 
Jupiters orbiting > 1.5 M stars . All radial velocity-d et ected planets around > 1.5 M ffi stars 
have semimajor axes > 0.5 AU Jjohnson et aI.ll2007al lbL liooS Isato et alll2008l : IWright et all 



2009|). 



Planetary migration may explain aspects of the semimajor axis distribution of Jovian 
planets, including the origin of hot Jupiters (which cannot form in situ ) and the period valley 
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for solar-mass stars (e.g. Ilda and Linl 12004 ; iBurkert and Idal 120071 ). However, the cause for 



the lack of hot Jupiters and other close-in giant planets around high-mass stars is less clear. 
One possibility is that hot J upiters surroundi ng high-mass stars are engulfed as the stars 
evolved off the main sequence (jSato et al.ll2008l ). However, most high-mass stars w ith planets 
(subgiants) are physically too small to engulf hot Jupiters (j Johnson et al.ll2007bl ). 

In this paper, I show that models of planet formation/migration with a stellar- mass- 
dependent lifetime of the gaseous circumstellar disks from which planets form can explain 
the dearth of hot Jupiters around high-mass stars as suggested by radial velocity surveys. 
If gas disks disappear much faster around high-mass stars than around solar and subsolar- 
mass stars, then inward migration is halted and the plan ets are stranded at large semimajor 
axes. The arguments described here build on work by IBurkert and Idal (120071 ) who stud- 
ied how a stellar-mass-dependent gas disk lifetime may explain why the period valley is 
more pronounced for planets orbiting F stars than G and K stars. In Section 2, I make 
simple analytical arguments to show that the observed semimajor axis distribution of giant 
planets may emerge from a stellar- mass-depend ent disk lifet i me. I n Section 3, I perform 



numerical modeling similar to recent work (e.g. Ilda and Linl l2004j ; IBurkert and Idal 120071 ) 
to test my hypothesis. I construct synthetic distributions of giant planets using a suite of 
Monte Carlo simulations with a range of stellar mass- dependent disk lifetimes and as well as 
a mass-independent lifetime. A mass-independent lifetime poorly reproduces the observed 
semimajor axis distribution of planets, while a mass-dependent lifetime reproduces observed 
trends. 



2. Analytical Motivation 



Planets form in disks around young stars. Once planets gr ow to Jovian masses , they 



can open a gap in the disk and undergo 'Type II migration' (ILin and Papaloizoul 11985 



Ward 



19971 ) with a migration rate regulated by the local viscous diffusion time. For regions 
interior to ~ 20 AU, migration is inward. The migration rate can be parameter i zed as suming 



a Minimum Mass Solar Nebula (MMSN) model as prescribed in Ilda and Linl (120041 ) scaled 
to the star's mass where the initial gas column density is £ s ,im ,mmsat = 2400 g cm -2 at 1 
AU: 



da 



— ~ 1.3 x 10~ 5 AU/yr x 
at 



a. 



\0.5/ 



J g,lMQ,MMSN 

M, 



1AU' 



M n ' v 10- 



' Mr. 



1.5 e -4t/r 9 _ 



(1) 



In this equation, M p /Mj is the planet mass in Jovian masses, a is the viscosity parameter, 
M* is the stellar mass, t is time, and r g is the timescale for the local disk surface density to 
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drop to < 1-2% of its original valuql 



The maximum Type II migration rate is the radial velocity of the gas, da p /dt ~ - 
1.5a(-^-) 2 a p Q, where H is the disk scale height at a p and Q is the Keplerian frequency For 
a MMSN surface density profile this rate is 



^~2.3x 10- s AUyr- 
dt 



\0.5i 



a 



10- 



(2) 



This equation governs the migration rate for disk masses that are much larger than the 
planet's mass. The nominal migration timescale from Equation 1, a p /(da p /dt), is then r m ji 
w 0.16 Myr for a Jovian mass planet at 1 AU with a = 5xl0~ 4 if r m> jj << r g . The 
timescale from the maximum drift rate (Equation 2) is about half this value. If the local 



surface density in a gaseous disk is drained to < 1% of its initial value (t 



), Equation 



1 implies a migration rate of < 1 AU/20Myr. Because nearly all gaseous circumstellar 
disks disappear by « 10 Myr f Currie et al. 2007al Jbh . migration is essentially halted. Thus, 
migration is halted if r g < r m ji. 

Simple arguments show that the gas disk lifetime may depend on stellar mass {r g oc 
M*~^) and that the ratio of r g to T m ji may strongly depend on stellar mass (see also 
Burkert and Idal 120071 ). Optical spectroscopic studies of 1-10 Myr-old star- forming regions 
reveals a stron g trend between accretion rate and stella r mass, dM*/dt oc M* 2 , for 0.03-3 
M Q -mass stars (jCalvet et al.ll2004j ; iMuzerolle et al.ll2005l ). The timescale for a star to accrete 
some fraction of its mass, x, is then r g ps 3Myr x (^j-)(-^) _1 or r g oc M* _1 . Thus, gas disks 
may dissipate faster around high-mass stars than around low-mass stars. 



A shorter gas dissipation timescale eventually leads to a slower migration rate. Ac- 
cording to Equation 1, at t/r 9i iM Q =0.2, a planet orbiting a 3 M star migrates inward at a 
slower rate than one orbiting a 1 M Q star if /3 w 1. At slightly later times (e.g., t/r s lMQ 
= 0.3), the migration rate for a planet orbiting 3 M star is half that for one orbiting a 1 
M star. Therefore, migration around higher-mass stars decelerates earlier; r g /r m ji will be 
smaller for higher-mass stars. Because smaller timescale ratios will strand more gas giants 
at larger semimajor axes (e.g. ~ 1 AU), the relative frequency of hot Jupiters and other 
close-in planets should be lower for higher-mass stars. 



1 I set the parameter f g in llda and Linl (|2004T ) equal to M^/M Q . Equation 65 in llda and Linl ([20041 ) is 



missing a factor of t/Mq/M*. 
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Numerical Model 



To test whether a stellar-mass-dependent gas disk lifetime can explain the observed 
distribution of extrasolar gas giant planets, I produce synthetic populations of exoplanets in 
semimajor axis vs. stellar mass s pace using a suit e of M onte Carlo simulations. My model 
generally follows the approach of iBurkert and Idal (120071 ) who show simulations for a 1 M Q 



star and vary the range in disk lifetimes (1-10 Myr, 3-30 Myr, 10-100 Myr) to show how a 
stellar mass- dependent disk lifetime regulates the distribution of exoplanets. Here, I perform 
simulations for a range of stellar masses and include an explicit power law dependence for 
the gas disk lifetime. 

Stellar masses for the synthetic population are randomized between 0.3 and 3 M Q with 
a probability distribution weighted towards solar /subsolar-mass stars (P(M*) oc M* -2,5 : a 
Salpeter-like IMF). To model the regions of planet formation, the planets' initial semimajor 
axes are chosen between 1 and 25 AU. For my fiducial model, I assume a gas disk lifetime 
that scales inversely with stellar mass (r g = 4 Myrx (M^/Mq) -1 ). 

For Jovian- mass planets to fo rm, I require that they reach an isolation core mass of 



Mj„ so = 5 M® (e.g. lAlibert et al.ll2005l ). I assume a 2.5 x scaled M MSN model fromllda and Lin 



(120041 ) to account for the median disk mass needed to form cores (IKenyon and Bromleyll2009l ) 
and set the metallicity comp arable to the media n metallicity of stars with detected Jovian- 
mass planets ([F e/H] ~ 0.15. IWright et all 120091) . I set the ice line location equal to values 
from Figure 1 of Kennedy and Kenyonl ( 2008 ) at 0.3 Myr interpolating between values for 
stars of different massed . Finally, I require that enough gas is left when the core mass is 5 
M e to form a Jovian- mass planet. This condition is equivalent to 



M 9iiso = IMMj x 



*0.75/ 



'2400#cm- 



(3) 



where the gas feeding zone size is 10 Hill radii. The Jovian planet-forming regions are > 2 
AU for 0.5 M stars, > 3 AU for 1 M stars, and > 5 AU for 2.5 M stars. 

From their birthplaces, I track the semimajor axis evolution of planets from Type II 
migration according to Equations 1 and 2, assuming a viscosity parameter of a=5xl0 -4 . 
Planets reaching within ~ 2-5 (~ 2-5 Rq(M*/M & ) ' 75 may be affected by magnetospheric 



disk truncation (ILin et al.lll996l ). which is not treated in this simple model. Any planet that 
reaches this small separation is given a final semimajor axis that is randomized between 2 



llda and Lin J 2004 ) and Burkert and Ida! (|2007l ) assume an optically-thin disk in determining the ice line 
position. Like Kennedy and Kenvon ( 2008h . I assume that the disk is not optically thin when planetesimals 



grow into gas giant cores. 
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and 5 R*. I make two simplifying assumptions in the fiducial model that will be removed 
later. First, 1 assume that all planets are the mass of Jupiter. Second, I assume that all 
planets form at 1 Myr: a characteristic time for the formation of 5 M cores at 5 AU around 



solar- mass stars (IKenyon and Bromley! 120091 ) 



Table 1 lists simulation results for a total of 20,000 planets (20 simulations of 1000 
planets; ^ total)- The table shows the total number of planets with final semimajor axes < 3 
AU and in three semimajor axis bins (< 0.2 AU, 0.2-0.5 AU, and 0.5-3 AU) for each of the 
stellar mass bins (0.3-0.5 M Q , 0.8-1.5 M , and 1.5-3 M ). The lower left panel of Figured] 
shows the final semimajor axis versus stellar mass distribution of these planets (black circles). 
I overplot the distribution of radial velocity-detected planets with well-constrained stellar 
masses (258; red star^j). The lower right panel shows a histogram plot of the semimajor 
axes for the synthetic population in three mass bins: 0.3-0.5 M (dashed line), 0.8-1.5 M Q 
(dotted line) and 1.5-3 M (solid line). 

Despite its simplicity, the model yields good agreement with observed semimajor axis 
distributions. In the synthetic population, stars with masses between 0.3 M and ~ 1 M 
have many hot Jupiters. The right panel of Figure [TJ shows that few 0.8-1.5 M stars 
have planets with a p ~ 0.1 and 0.5 AU, consistent with the 'period valley' in the observed 
population. 

Most strikingly, the synthetic population shows a sharp drop in the number of hot 
Jupiters from ~ 1.3 M to 1.5 M . For my model assumptions, the synthetic population lacks 
any hot Jupiters for stars with M* > 1.6 M and lacks any planets with intermediate (0.2-0.5 
AU) distances for M* > 1.8 M . This trend is consistent with the clear lack of detected 
hot Jupiters around 1.5-3 M stars. The x 2 values comparing the relative frequencies of 
hot Jupiters, planets at 0.2-0.5 AU, and those at 0.5-3 AU with the observed distribution 
confirm that agreement is good, especially for > 0.8 M stars. 

The semimajor axis distribution shows far poorer agreement if the gas disk lifetime is 
independent of or weakly dependent on stellar mass (Figure [2]). Assuming r g =4 Myr for all 
masses turns all planets around > 1.5 M stars into hot Jupiters, weakens the period valley 
for solar-mass stars, and confines all planets around low-mass stars to semimajor axes with a p 
> 10 AU (top-left panel). These properties are inconsistent with the observed distribution. 
The model with r g = 2 Myr (top-right panel) for all stars yields the correct distribution 
for high- mass stars but eliminates all hot Jupiters around solar/subsolar-mass stars, and 
strands all planets around subsolar-mass stars at a p > 10 AU. Models with r g oc M* -0,5 fare 
marginally better (bottom panels). The model with t 9) im = 4 Myr predicts a pileup of hot 
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Jupiters, many planets at ~ 3-10 AU, and a dearth of planets at ~ 0.2-0.5 AU for all stars 
with masses > 0.8 M . It also confines planets around < 0.5 M Q stars to > 0.3 AU with a 
peak at ~ 10 AU. The model with T g) iM e = 2 Myr confines all planets to a p > 3 AU. These 
features are clearly not present in the observed population. The x 2 values for all of these 
simulation runs exceed 100 for low-mass stars and (sometimes) high-mass stars. 

Motivated by the success of the fiducial model with r g oc M* -1 , I remove assumptions 
regarding the planet's mass and core formation timescale and run a second set of simula- 
tions. First, I set the planet's mass equal to the minimum gap-opening mass for Type II 
migr ation, requiring that the planet's Hill radius is larger than the local disk scale height 
(e.g. 



Lin and Papaloizoul 1985 ): 



Second, I estimate the formation timescale for each core explicitly, require that the core 
can form before gas is dissipated, an d require that enough gas is left to form a migrating 



planet after core formation. From the iKenyon and Bromleyl (120091 ) results, I extrapolate the 
formation timescale of 1 Myr at 5 AU around a solar-mass star with = 2.5 g cm -2 to 
different distances, stellar masses, and dust column densities. I assume that r mrf , oc (S^)" 1 , 



which is a reasonable approximation of the numerical results ( IKenyon and Bromleyl 120081 ) 



Results from the second set of simulations (the bottom half of Table 1, Figure ED show 
good agreement between the synthetic distribution and the observed distribution if r g oc 
M* -/3 , where (3 = 0.75-1.5. The semimajor axis versus stellar mass distribution (shown 
for T g =2 Myrx M* _1 , top-left panel) features the same sharp drop in the frequency of hot 
Jupiters and planets with a p = 0.2-0.5 AU (log(a p ) = -0.7- -0.3) for M* > 1.5 M exhibited 
by the fiducial model. The histogram plots for (3=1 (top right panel) and (3 = 0.75 (second 
row, left panel) confirm that there is a lack of hot Jupiters and other planets with a p < 
0.5 AU for high- mass stars. The models successfully reproduce the two other features of 
the observed distribution: the period valley for solar-mass stars and the lack of long-period 
planets around low-mass stars. The agreement is confirmed quantitatively as models with (3 
= 0.75-1.5 typically have \ 2 values less than ~ 20 for all stellar mass bins. 

I also calculate the percentage of planets that achieve a gap-opening mass and undergo 
Type II migration (N 9ap ) for each of the stellar mass bins (100xN 9ap /N to t a ;). For (3 < 1.5, the 
frequency is lowest for M stars. This percentage probes the relative sizes of planet-forming 
regions for a given scaled dis k mass (2.5 x a scaled MMSN) and is different from but related 



to the frequency of planets (I Johnson et al.ll2007bl ). For calculations with lower scaled disk 



masses (e.g., 1.25 x MMSN scaled) the planet-forming regions shrink, and only solar-to-high 
mass stars form migrating planets. For calculations with an even lower scaled disk mass, 
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only high-mass stars form migrating planets. The percentage of migr ating gas giants for M 
stars is lowest because their disks have lower masses: M^fc oc f g oc M* (Kennedy and Kenyon 



20081 ). 



Models with r G = 2 Myr and (3 < 0.5 (third row) begin to show disagreement with the 
observed distribution as they predict a pileup of hot Jupiters and period valley for high- 
mass stars and a lack of hot Jupiters for low-mass stars. The x 2 values for low-mass stars 
and high-mass stars for (3 < 0.5 are significantly higher than for (3 = 0.75-1.5. Therefore, 
according to my simulations, the observed semimajor axis distribution of gas giant planets 
results from a stellar-mass-dependent gas disk lifetime. 



4. Discussion 



Through a series of Monte Carlo simulations, I have shown that a stellar-mass-dependent 
gas disk lifetime can explain the observed semimajor axis distribution of extrasolar gas giant 
planets, including the lack of hot Jupiters and other planets with a p < 0.5 AU around high- 
mass stars. Synthetic distributions of planets produced assuming that the gas disk lifetime, 
r g , scales as M* _/3 with (3 = 0.75-1.5 reproduce key features in the observed semimajor 
axis distribution. Distributions from models lacking a stellar-mass-dependent disk lifetime 
quantitatively provide a poorer match to observations. 



This work extends and complements the investigation of iBurkert and Idal (120071 ) who 
use Monte Carlo simulations and semi-analytical prescriptions for planet growth to explain 
exoplanet tre nds for solar-mas s stars . While high-mass stars lack a pronounced period valley, 
predicted by IBurkert and Idal (120071 ). the dearth of planets at 0.2-0.5 AU agrees with their 
pr edictions. Moreoy er, this work shows that a stellar mass-dependent disk lifetime, invoked 
by IBurkert and Idal to explain exoplanet trends for solar-mass stars, may explain trends for 
planets around stars with a wide range of masses. Future modeling work is necessary to 
test this hypothesis more conclusively. Future modifications include incorporating a more 
sophisticated treatment of circumstellar gas accretion, modeling Type I migration, deter- 
mining the sensitivity of the planets' synthetic distributions to the migration r ate (i.e., value 



of a) , and tracking the migration of planets while they are accreting gas (e.g. lAlibert et al. 



20051 ). 



Rec ent studies of young stars i n clusters support a stellar mass-dependent gas disk 
lifetime (Kennedy and Kenyonl 120091 ) . Based on optical spectroscopy, the frequency of gas 
accretion in 2-15 Myr-old clusters is significantly higher fo r stars with M* < 1 M^ than for 
higher-mass stars (e.g., IC 348, Tr37, and h and \ Persei iDahml 120081 ; ICurrie and Kenyon 
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20091 : ISicilia-Aguilar et alJliooS ICurrie et al. l l2007bh . Secondary characteristics of gas-rich 
disks (optically thick thermal infrared em ission) are also rarer for high-mass stars (e.g., 
Carpenter et al.ll2006l ; ICurrie et al. Il2009bl ). Combining cluster data to empirically constrain 
T g (M ic ) may be possible. However, uncertainties in stellar ages for stars in the youngest (< 
3 Myr) clusters present a strong challenge to constructing an empirically based gas disk 
lifetime. I will address these issues in a future paper. 

If the gas disk lifetime strongly depends on stellar mass, my model simulations suggest 
that future radial velocity surveys will find few gas giant planets orbiting at small separations 
from high-mass stars. If ~ 10% of high-mass stars have gas giant planets with a p < 3 AU 
( I Johnson et al.ll2007bl ). the simulations with f3 = 0.75-1.5 imply that out of 1,000 high-mass 
stars targeted for radial velocity surveys, fewer than w 15 will have planets at a p < 0.5 AU 
while more than ss 75 will have planets at a p > 0.5 AU. Ongoing surveys will provide a larger 
sample from which to compare observed and predicted frequencies of hot Jupiters and other 
gas giants at small separations (J. Johnson, in preparation). 



I thank John A. Johnson and Jason Wright for discussions on exoplanet data and for 
encouraging the writing of this paper; Geoff Marcy, Scott Kenyon, Dan Fabrycky, Ruth 
Murray-Clay, Brad Hansen, and Jonathan Irwin for valuable comments; and the referee for 
a rapid and insightful report. This work is supported by NASA Astrophysics Theory Grant 
NAG5-13278 and NASA Grant NNG06GH25G. 
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Table 1. Simulation Results 



ID t p Number Fractions of Planets in Each Semimajor Axis Bin \ 2 % 

a p <0.2 AU a p =0.2-0.5 AU a p =0.5-3 AU 
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IM 


HM 


LM 


IM 


HM 


LM 


LM 


HM 


LM 


IM 


HM 


LM 


IM 


HM 


LM 


IM 


HM 


Obs. 








11 


176 


20 


0.73 


0.27 


0.0 





18 


0.13 


0.0 


0.09 





60 


1.00 








1.8" 


4.2° 


8.9 a 


Ml 


4 


1 




5524 


2198 925 


0.50 


0.42 


0.05 





09 


0.10 


0.05 


0.41 





49 


0.89 


28.8 


6.3 


3.0 










4 





5 


1957 


2469 3499 


0.10 


0.48 


0.52 





11 


0.09 


0.09 


0.79 





43 


0.39 


157.0 


13.1 


114.5 










4 







602 


2787 4782 





0.52 


0.94 







0.09 


0.03 


0.99 





39 


0.04 


245.2 


19.0 


317.9 










2 





5 


329 


255 


20 






















1.00 


1 


00 


1.00 


245.2 


44.0 













2 







92 


365 


1290 








0.14 










0.09 


1.00 


1 


00 


0.77 


245.2 


44.0 


14.2 








M2 


2 


1 




948 


893 


1146 


0.64 


0.46 


0.06 





05 


0.10 


0.10 


0.32 





11 


0.84 


13.7 


11.0 


6.9 


13.4 


20.8 


20.9 




2 


1 


5 


1337 


629 


234 


0.76 


0.48 


0.0 





04 


0.10 


0.03 


0.2 





41 


0.97 


5.7 


14.3 


0.3 


19.6 


21.4 


14.5 




2 





75 


783 


911 


1222 


0.53 


0.52 


0.15 





07 


0.10 


0.10 


0.40 





38 


0.75 


25.1 


19.7 


16.7 


11.1 


20.7 


22.4 




2 





5 


546 


584 


780 


0.40 


0.55 


0.24 





08 


0.09 


0.13 


0.52 





36 


0.63 


53.5 


24.2 


37.2 


9.2 


20.4 


26.1 




2 





1 


319 


671 


975 


0.31 


0.53 


0.51 





09 


0.08 


0.12 


0.60 





40 


0.36 


78.3 


19.4 


120.4 


5.8 


22.0 


30.3 




3 


1 




1173 


810 


751 


0.74 


0.65 


0.26 





05 


0.08 


0.12 


0.21 





27 


0.62 


5.5 


45.0 


39.9 


18.1 


26.6 


25.8 




4 


1 




1499 


961 


998 


0.78 


0.78 


0.46 





04 


0.04 


0.10 


0.18 





18 


0.45 


5.3 


78.3 


92.2 


21.8 


29.2 


30.6 



Note. — Statistics listing the distribution of low- mass stars (LM; 0.3-0.5 Mg), intermediate-mass stars (IM; 
0.8-1.5 Mg), and high-mass stars (HM; 1.5-3 Mg) for the observed population (Obs.) and the two model 
runs (Ml and M2). "Number" refers to the number of planets with final semimajor axes < 3 AU. The relative 
fraction of planets is computed for a p <0.2 AU, a p =0.2-0.5 AU, and a p =0.5-3 AU for each stellar mass bin. The 
X 2 parameter is computed by comparing the relative fraction of simulated planets with the observed fraction for 
each semimajor axis bin for a sample size equal to the number of observed planets in each stellar mass bin. The 
column "%" lists the percentage of gap-opening planets for each stellar mass bin regardless of final semimajor 
axis. Percentages are not calculated for the first set of simulations because in those simulations planets are 
assumed to form whenever there is sufficient mass for them to form, (a) — Shown for reference is the absolute 
frequency of all detected planets from radial velocity surveys jjohnson et al.ll2007bT) for different stellar mass 
ranges. 
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-2_0 -1.5 -10 -0.5 0.0 0.5 1.0 1.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

log(a/1AU) log(a/1AU) 



Fig. 1. — Top panels: the observed semimajor axis vs. stellar mass distribution (left) and 
the histogram plot of semimajor axes from radial velocity surveys. Bottom panels: the 
semimajor axis vs. stellar mass distribution for the synthetic population (black circles) with 
my fiducial model, assuming r g oc M* -1 (left) and distribution of radial velocity-detected 
planets (red stars), (right) Histogram plots of the semimajor axis distribution for 0.3-0.5 
M stars (dashed line), 0.8-1.5 M stars (dotted line), and 1.5-3 M stars (solid line). 
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Fig. 2. — Histogram plots of the semimajor axis distribution assuming r g = 4Myr (top 
left), r g = 2 Myr (top right), r g = 4 Myrx (M*/M )-°- 5 (bottom left), and r g = 2 
Myrx (M*/M ) _a5 (bottom right) for my fiducial model. Each distribution is inconsistent 
with the observed population. 
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Fig. 3. — Results for the second set of simulations. Top left: distribution of semimajor axis 
vs. stellar mass and histogram plot of the semimajor axis distribution (top right) for r g = 
2 Myrx (M^/Mq)" 1 . Second row: histogram plots for r g = 2 Myrx (M*/M ) a75 (left) and 
r g =3 Myrx(M^/M )~ 1 (right). Bottom row: histogram plots for r g = 2 Myrx (M*/M )-°- 5 
(left) and r g = 2 Myrx (M*/M )- ai (right). 



